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SUMMARY 


This is the final report for contract DAAKO2-70-C-0258 and covers 
the period from January 1971 through June 1971. !t is a summary of ali 
programs written and/or modified under this contract during that period by 
Braddock, Dunn and McDonald, Inc. for the U.S. Army Mobility Equipment 
Research and Development Center and subsequently, for U.S. Army Harry 
Diamond Laberatories. Each new program is described in the form of a user's 
guide which includes a discussion on how the program operates, instructions on 
how to execute the program, and sample output generated by the program. 


For each modified program the discussion tells how and why it was modified. 


FOREWORD 


The programs documented in this report are part cf a comprehensive 
package of digital computer codes which calculate the electromagnetic 
pulse (EMP) environment that is produced by a nucivar explosion. This 
work is a continuation of an extensive technical effort formerly conducted 
at the U.S. Army Mobility Equipment Research and Development Cenier 
(MERDC), Fort Belvoir, Virginia, and presently conducted at the U.S. Army 
Harry Diamond Laboratories, Washington, D.C. The modification of previously 
documented computer codes and the writing of new codes has necessitated this 
documentation. 

The programming tasks were performed under contract DAAKO2-70-C-0258 
by Mr. Jeffrey A. Borbely (project leader), Mr. David L. Jones, and 
Mrs. Deborah H. Stump of Braddock, Dunn and McDonald, Inc. Technical 
supervisor was Mr. William T. Wyatt, Jr. formerly of the Physics Division, 


Electromagnetic Effects Laboratory, MERDC, and presently of Branch 1030, HDL. 
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|. NEW CODES AND MAJOR MODIFICATIONS 
A. COMPARE 
1. Introduction 
Program COMPARE is a plot program which draws from one to ten 
parallel curves on the same plotting area to facilitate data comparison. 
Given a set of EMP field values which are a function of three independent 
parameters (range, time, and theta), the user may elect to hold any one of 
the three parameters constant, have each of the ten curves represent a 
different value of the second parameter, and plot the EMP field values 
versus the remaining parameter. Each of the curves may be scaled inde- 
pendently or the data for al] the curves may be scaled as one set. 
2. Plot Definition and Control 
The EMP field data is supplied to program COMPARE by reading the 
output tapes from program ELECTRA the ground burst EMP cogs. 8 In arder 
to determine which data is to be read from the tape a series of { iput vata 


cards are read by program COMPARE. The detailed structure of these cards 


1. Jones, D. L. and D. H. Stump, ELECTRA, ORESTES, and Supporting Graphic 
Display Codes, U.S. Army Mobility Equipment Research and Development 
Center, Fort Belvoir, Virginia, February 1971, pp. 1-36. 


2. Borbely, J.A. and D. L. Jones, ELECTRA, An Electromagnetic Pulse 
Fortran Program (User's Guide), U.S. Army Mobility Equipment Research 
and Development Center, Fort Belvoir, Virginia, October 1969. 
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will be discussed in Section 1.A.10,. The first data card determines which 
type of plot is desired. That is, it establishes which parameter wil! be 
hetd constant, which parameter (designated z-axis’ will change only from 
curve to curve, and which parameter (designated x-axis) will be plotted 
against the field data. This first card will contain one of the following 


ten-character keywords: 


PLOT TYPE KEY WORD Z-AXIS X-AXIS CONSTANT 
] TIME-THETA TIME THETA RANGE 
2 RANGE-THETA RANGE THETA TIME 
3 Ti ME-RANGE TIME RANGE THET#. 
4 RANGE-TIME RANGE TIME THETA 
5 THETA-T tME THETA TIME RANGE 


These five key words are used as switches to control the execution cf 
five separate sections of program COMPARE. The THETA-RANGE-TIME combina~ 
tion is not implemented in this version of COMPARE. 

After the program enters one of the five sections, the additional 
data items needed to set up that type of plot are read from cards. These 
data items include the value of the constant parameter, the values of the 
2 axis parameters, and the values or limits of the x-axis parameters, In 
additton controls for scaling and for x-axis spacing are also set by data 
card. Once these parameters have been set, the EMP data is read from tape, 
edited and ordered accordingly, and then transferred to the appropriate 


plotting subroutine (LINEUP, LYNEUP, LTCONE, LTKONE, TIMAXS, or TYMAXS) . 


3. Subroutines LINEUP and LYNEUP 

Subroutine LINEUP is called by program COMPARE to produce plots 
of types 1 and 2 as defined in Section 1.4.2. It calls subroutine SCALE, 
treating the EMP field data for all curves as one set, and receives one 
maximum, one minimum, and one delta y scale factor. LINEUP then draws each 
curve separately, scaling each data set against the maximum, the minimum, 
and the delta y scale factor. Each curve has its own set of axes 
appropriately annotated. 

If the plot is of type 1, the constant parameter is range and is 
annotated by the statement ALL CURVES ARE FOR R = (range). Each curve 
reprvsents a different value of either source time or retarded time and 
is labeled either T = (source time) or Tt = (retarded time). Tne x-axis 
shows the values of theta which may be spaced according to actual value 
or spaced at equal intervals as a user option. 

If the plot is of type 2, the constant parameter is either source time 
or retarded time and is annotated either by the statement ALL CURVES ARE 
FOR T = (source time) or by ALL CURVES ARE FOR t = (retarded time). Each 
curve represents a different value of range and is labeled R = (range). 
The x-axis again shows the values of theta as described above. 

Subroutine LYNEUP differs from LINEUP only in the area of data 
scaling. In subroutine LYNEUP the EMP field data set for each curve is 
supplied to subroutine SCALE separately. The resulting maximum, minimum, 
and scale factor fer each curve are then reflected in the annotation on its 


set of axes, 


4. Subroutines LTCONE and LTKONE 
Subroutine LTCONE is called by program COMPARE to produce plots 
of type 3 as defined in Section 1.A. 2. Im a manner similar to that 
described for subroutine LINEUP in Section 1.A.3., subroutine LTCONE 
calls subroutine SCALE, treating the EMP field data for all curves as one 
set, and then draws the curves separately with individually annotated 
axes. 

For a type 3 plot the constant parameter is theta and is annotated 
by the statement ALL CURVES ARE FOR 9 = (theta). Each curve represents 
a different value of source time ard is labeled T = (source time). The 
x-axis shows values of retarded range R’ where R'=R-cxT, R = actual 
range, T -- source time, and c = speed of light. Since each curve 
represents a different value of T, different retarded ranges must be 
calculated for each curve and then each curve must be shifted to align 
with the annotated x-axis at the bottom of the plot. 

Subroutine LTKONE differs from LTCONE only in the area of data 
scaling. In subroutine LTKONE the EMP field data set for each curve is 
supplied to subroutine SCALE separately. The resulting maximum, minimum, 
and scale factor for each curve are then reflected in the annotation on 
its set of axes. 

5. Subroutines TIMAXS and TYMAXS 

Subroutines TIMAXS and TYMAXS are called by program COMPARE to 


produce plots of types 4 and 5 as defined in Section 1.A.2. Subroutine 


TEMAXS calls subroutine SCALE, treating the EMP field data for al! curves 

as one set, and receives one maximum, one minimum, and one delta y scaling 
factor. TIMAXS then draws each curve separately, scaling the data against 
the maximum, minimum, and scale factor established for the total data set. 
Each curve has its own set of axes appropriately annotated. 

Subroutine TYMAXS supplies subroutine SCALE with the data for each 
curve separately and then plots each curve according to its own maximum, 
minimum, and deita y scale factor. Each curve again has its own set of 
axes appropriately annotated. 

If the plot is of type 4, the constant parameter is theta and is 
annotated by the statement ALL CURVES ARE FOR 6 = (theta). Each curve 
represents a different value of range and is labeled R = (range). The 
x-axis shows the values of either source time or retarded ‘ime in units of 
nanoseconds and is labeled either T IN NSEC or t IN NSEC. The time values 
may be spaced according to actual value or at equal intervals as a user 
option. 

If the plot is of type 5, the constant parameter is range and is 
annotated by the statement ALL CURVES ARE FOR R = (range). Each curve 
represents a different value of theta and is labeled 6 = (theta). The 
x-axis is the same as described above for type 4 plots. 

6. Subroutine COEF 

Subroutine COEF helps set up the y-axis labeling for the six 


plotting subroutines described in Sections 1.A.3., |1.A.4., and 1.A.5. 


Snteead 


Given a delta y scale factor, subroutine COEF breaks the factor apart into 
a coefficient between one and ten and a power of ten. This coefficient is 
then used to determine the number of intervals into which the y-axis will 
be divided. The power of ten is used along with the minimum y value to 
determine the annotation to be written beside each tic mark on the y-axis. 
7. Subroutine SCALE 
Subroutine SCALE performs data scaling for the six plotting 
subroutines described in Sections 1.A.3., !'.A.4., and 1.A.5. Given a data 
array, subroutine SCALE scans the array to find the maximum and minimum 
values. {t then defines the delta y scale factor by the following 
expression: 
Ay = (maximum y - minimum y)/ length of plot axis 
8. Subroutine LINES 
Subroutine LINES is a plot routine designed for drawing data 
curves. Given two data arrays which are to be plotted against each 
other, subroutine LINES performs a curve fit to the data in a manner such 
that the curve can be represented by a series of equally spaced dots, line 
segments or both. 
9. Subroutines GREEK, SEPAR8, and {SHIFT 
These three subroutines work togcther to draw the letters of the 
Greek alphabet in either upper case or iower case form. They wil! be 
described in detail in Section |.C. of this report. They are called 


by program COMPARE to supply t's and 8's for labeling the plots. 
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10. Description of tnput 
The card input for program COMPARE can best be described by 
breaking it up into five distinct groups ev2n though several of the 
groups will hav2 the same card types. 
a. Type 1 
The input for plots of type | as defined in Section 1.A.2. 


consists of the five card types (1A, 1B, IC, 1D, and IE) shown below: 


TYPE DESCRIPTION 


Ea CARD VARIABLE NAME 
Ea COL. OR CONTENTS 


| PROGRAM 
NAME 


1A Select type of CHOICE |Alpha-|Must equal TIME-THETA 
plot betic 
11-80] Not used Blank 
Number of NZ Inte- {Determines number of 
curves ger time values to appear 
on card IC below. Must 
be < 10 in present 
version 
11-80] Not used Blank 
1c Time value Z(1) Float-|May be source time or 
for Ist curve ing retarded time in seconds 
(See TCOL below) 
11-20] Time value Z(2) Float-|May be source time or 
for 2nd curve ing retarded time in seconds 
(See TCOL below) 
21-30| Time value Z(3) Float-|May be source time or 
for 3rd curve ing retarded time in seconds 
(See TCOL below) 
31-40] Time value Z(4) Float-|May be source time or 
for 4th curve ing retarded time in seconds 
(See TCOL below) 


for 5th curve ing retarded time in seconds 


41-50| Time value Z(5) Float-|May be source time or 
(See TCOL below) | 


CARD | CARD VARIABLE NAME 
TYPE | COL. OR CCNTENTS 


Time value 
for 6th curve 


Time value 
for 7th curve 


Time value 
for 8th curve 


Time 
increment 


Range 
increment 


Constant 
parameter 


Time 
indicator 


Length of 
x-axis 


Length of 
y-axis 


Not used 


EMP field 
to be plotted 


Not used 


Minimum 6 
value 


PROGRAM 
NAME 


CONVAR 


TCOL 


XLEN 


YLEN 


THETMN 


TYPE 


DESCRIPTION 


May be source time or 
retarded time in seconds 
(See TCOL below) 


May be source time or 
retarded time in seconds 
(See TCOL below) 


May be source time or 
retarded time in seconds 
(See TCOL below) 


Used as tolerance factor 
in matching time values 


Used as tolerance factor 
in matching range values 


Range value for all 
curves 


TCOL=0, retarded time 
TCOL40, source time 


6 inches recommended 


This is the combined 
total length of all the 
y-axes. Each individual 
axis length will equal: 
(YLEN/NZ) - .25 

YLEN must be < 10. 


Blank 


Must equal BPHI, ETHE, 
ERAD, or SIGT. 


Blank 


Minimum theta value for 
x-axis data. 
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CARD | CARD DESCRIPTION 


VARIABLE NAME | PROGRAM | TYPE 
QR CONTENTS NAME | 


THETMX Float- |Maaimum theta value for 
ing  |x-axis data. 


Maximum 6 
value 


[tCX Inte- 'INCS=1, spaced accord- 
ger ling to data; INCX=2, 


equally spaced 


x-axis spacing 
control 


y-axis scaling {SCALE |Inte- 
control | ger 


| 


curves as one Set; 
ISCALE #40, scale each 
curve separately. 


Blank 


| 

{ 

| 

ISCALE=0, scale all 

Not used | 


A standard deck for type | plots will consist of one type IA, one type 1B, 
one type 1C (a second type IC is required if NZ is greater than eight), 
onc type 1D, and as many type IE cards as desired. The last type JE card 
must be followed by a blank card which signals the end of that set of type 
1 plots. 
b. Type 2 
The input for plots of type 2 as defined in Section |.A.2., 


consists of the five card types (2A, 2B, 2C, 2D, 2E) shown below: 


CARD | CARD VARIABLE NAME PROGRAM | TYPE DESCRIPTION 
TYPE | COL. OR CONTENTS 


Select type Must equal RANG-THETA 
of plot 


Not used Blank 


Number of Determines number of 
curves range values to appear 


on card 2C below. Must 
be < 10 in present version, 


re TOTO OCU TOOPO 


DESCRIPTION 


CARD } CARD VARIABLE NAME PROGRAM | TYPE 
TYPE | COL. OR CONTENTS 


26 11-80 | Not used Blank 


cont. 
2c 1-10 | Range value Range in meters 
for Ist curve 
11-20 | Range value Range in meters 
for 2nd curve 
21-36 | Range value Range in meters 
for 3-d curve 
31-40 | Range value Float- {Range in meters 
for 4th curve 
41-50 | Range value Range in meters 
for 5th curve 
51-60 | Range value Range in meters 
for 6th curve 
61-70 | Range value Range in meters 
for 7th curve 
71-80 | Range value Range in meters 
for 8th curve 
2D 1-10 | Time DT Float- [Used as tolerance factor 
increment ing in matching time values 
11-20 | Range DR Float- |Used as tolerance factor 
increment ing in marching range values 
21-30 | Constant CONVAR |Float-|Time value for all curves, 
parameter ing May be source time or 
retarded time 
31-40 | Time TCOL Float~ ITCOL=0, retarded time 
indicator ing TCOL#0, source time 
41-50 | Length of XLEN Float- |6 inches recommended 
x-axis | 
10 


CARD | VARIABLE NAME 


COL. ; OR CONTENTS 


DESCRIPTION 


PROGRAM | TYPE 
NAME 


51-60 | Length of YLEN 


y-axis 


Float-|This is the combined 
ing total length cf all the 
y-axes. Each individ- 
ual axis length will 
equal: (YLEN/NZ) -.25 
YLEN must be < 10. 


Not used 


EMP field 
to be plotted 


Blank 


Must equal BPHI, ETHE, 
ERAD, or SIGT. 


Not used Blank 


Minimum 6 THETMN 


value 


Minimum theta value for 
x-axis data. 


Maximum 6 
value 


THETMX Maximum theta value for 


x-axis data. 


x-axis spacing 
control 


INCX INCX=1, spaced according 
to data; INCX=2, 

equally spaced 

y-axis scal- iSCALE 


ing contro] 


ISCALE=0, scale all 
curves as one set; 
ISCALE # 0, scale 

each curve separately. 


Blank 


A standard deck for type 2 plots will consist of one type 2A, one type 

2B, one type 2C (a second type 2C is required if NZ is greater than eight), 
one type 2D, and as many type 2E cards as desired. The lest type 2E card 
must be followed by a blank card which signals the end of that set of type 


2 plots. 


c. Type 3 
The input for plots of type 3 as defined in Section 1.A.2. 


consists of the six card types (3A, 3B, 3C, 3D, 3E, and 3F) shown below: 
CARD | CARD 
TYPE | COL. 


ae 
11-80 


VARIABLE NAME 
OR CONTENTS 


PROGRAM 
NAME 


DESCRIPTION 


CHOICE 


Select type 
of plot 


Must equal T!ME~-RANGE 


Not used Blank 


Determines number of 
time values to appear 

on card 3C below. Must 
be < 10 in present version. 


Number of 
curves 


Not used Blank 


Time value z(1) Float-|Source time value in 
for Ist curve ing seconds 

11-20 | Time value Z(2) Float-|Source time value in 
for 2nd curve ing seconds 

21-30] Time value Z(3) Float-|Source time value in 
for 3rd curve ing seconds 

31-40] Time value Z(4) Float~ {Source time value in 
for 4th curve ing seconds 

41-50] Time value Z(5) Float~|Source time value in 
for 5th curve ing seconds 

51-60 | Time value 2 (6) Float-|Source time value in 
for 6th curve ing seconds 

61-70] Time value Z(7) Float-|S<  e time value in 
for 7th curve ing secunds 

71-80 | Time value z(8) Float-|Source time value in 
for 8th curve ing seconds 


Seen a Ree ee nl ee 


3E 


3F 


CARD 
COL. 


21-36 


31-40 


Time 
increment 


Range 
increment 


Length of 
x-axis 


Length of 
yraxis 


Not used 
Number of 


range values 


First range 
value 


Increment 
between range 
values 

Not used 


EMP field to 
be plotted 


Not used 


Cons tant 
parameter 


X-axis spacing 
control 


VARIABLE NAME PROGRAM | TYPE DESCRIPTION 
OR CONTENTS NAME 


DT Used as tolerance factor 
in matching time values 
DR Used as tolerance factor 
in matching range values 
XLEN 6 inches recommended 
YLEN This is the combined 


total length of al! the 
y-axes. Each individual 
axis length will equal 
(YLEN/NZ) - .25 YLEN 
must be < 10, 


Blank 


NRANGE |Inte- |Number of range values 


ger to be calculated for 
X-axis. 
RMIN Float-|First range value, basis 
ing for calculating all other 


range values 
RINC Float- |Each range value wil] 
ing equal the previcus 
| value plus RINC 
Blank 


AVAR Alpha- jMust equal BPHI, ETHE, 
betic {ERAD, or S!GT 


Blank 


CONVAR Float- |Theta value for all 


ing curves 
INCX Inte- |INCX=1, spaced accord- 
ger ing to data; INCX=2, 


equally spaced 


VARIABLE NAME 
OR CONTENTS 


PROGRAM | TYPE DESCRIPTION 
NAME 


ISCALE |tnte~ | ISCALE=0, scale all 
curves aS one set; 
ISCALE # 0, scale each 
curve separately. 


CARD ,; CARD 
TYPE | COL. 


control 


41-80 | Not used Blank 


A standard deck for type 3 plots will consist of one type 3A, one type 


38, one type 3C (a second type 3C is required if NZ is greater than eight), 
one type 3D, one type 3E, and as many type 3F cards as desired. The last 
type 3F card must be followed by a blank card which siqnals the end of that 
set of type 3 plots. 
d. Type 4 
The input for plots of type 4 as defined in Section {.A.2. 


consists of the seven card types (4A, 4B, 4C, 4D, 4E, 4F, and 4G) shown 


TYPE 


Alpha- 
betic 


below: 


| caro | carb | VARIABLE NAME 


TYPE | COL. OR CONTENTS 


PROukAM DESCRIPTION 


NAME 


Must equal RANGE-TIME 


Select type CHOICE 


of plot 


Blank 


Not used 


Determines number of 

range values to appear 
on card 4C below. Must 
be < 10 in present version, 


Number of curves 


Blank 


Not used 
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VARIABLE NAME 
OR CONTENTS 


PROGRAM 
NAME 


CARD | CARD DESCRIPTION 
TYPE | COL. 


1-10 


Float- 
i 


Range value 
for Ist curve 


Range i 


Range value Z(2) | Float- Range in meters 
| 


for 2nd curve ing 
for 3rd curve ing 


Z(4) 


Float- |Range in meters 
ing 


Range value 
for 4th curve 


Z(5) 


Float- |Range in meters 
ing 


Range value 
for 5th curve 


z (6) 


Float- |Range in meters 
ing 


Range value 
for 6th curve 


Float- 
ing 


Range value Z(7) Range in meters 


for 7th curve 


z(8) Float- 


ing 


Range value Range in meters 


for 8th curve 


Range value Z(3) Float \Ranoe in meters 


Used as tolerance factor 
in matching time values 


Time increment 


Used as tolerance factor 
in matching range values 


Range increment | DR 


Length of XLEN 6 inches recommended 


x-axis 


This is the combined 
total length of all the 
y-axes. Each individual 


Length of YLEN 
y-axis 


(YLEN/NZ) -.25 YLEN 
must be < 10. 


Blank, 


Not used 


CARD | CARD VARIABLE NAME PROGRAM ; TYPE DESCRIPTION 
TYPE | COL. OR CONTENTS NAME 


1-10 Number of NDELTS 
| different DELT's 


This controls the number 
of type 4F cards to be 
read. 


Minimem time 
value 


First time value, basis 
for calculating ell 
other time values. May 
be source or retarded 
time (See TCCL below). 


Not used Blank 


LF 1-10 | Number of times Each time value for the 
to use the DELT x-axis is equal to the 
value found on previous value plus 
this card DELT. NDELT tells how 

many times to add DELT. 
11-20] Increment See description of 
between time NDELT above. 


values. 


Not used Blank 


EMP field to AVAR Alpha-|Must equal BPHI, ETHE, 
be plotted betic |ERAD, or SIGT. 


Not used Blank 


Constant CONVAR |Float-| Theta value for all 
parameter ing curves 


x-axis spacing INCX Inte- | INCX=1, spaced according 
control ger to data; INCS=2, 
equally spaced. 


y-axis scaling ISCALE |Inte- | ISCALE=0, scale all 

control ger curves aS one set; 
ISCALE # 0, scale each 
curve separately. 


Time indicator TCOL Float-| TCOL=0, retarded time 
ing TCOL40, source time 


Not used Blank 


ot 


A standard deck for type 4 plots will consist of one type 4A, one type 
4B, one type 4C (a second type 4C is required if NZ is greater than eight), 
one type 4D, one type 4E, as many type 4F as specified by the NDELTS 
parameter, and as many type 4G cards as desired. The last type 4G card 
must be followed by a blank card which signals the end of that set of type 
4 plots. 
e. Type 5 
The input of plots of type 5 as defined in Section |.A.2. 


consists of the seven card types (5A, 5B, 5C, 5D, 5E, 5F, and 5G) shown 


be low: 
CARD | CARD VARIABLE NAME PROGRAM | TYPE | DESCRIPTION 
TYPE | COL. OR CONTENTS NAME 
5A 1-10 | Select type CHOICE | Alpha-|Must equal THETA-TIME 
of plot betic 
11-80] Not used Blank 
5B 1-10 | Time increment DT Float-|Used as tolerance factor 
ing in matching time values 
11-20} Range increment | DR Float-|Used as tolerance factor 
ing in matching range values 
21-30] Length of XLEN Flota-|6 inches recommended 
x-axis ing 
31-40] Length of YLEN Float-|This is the combined 
y~axis ing total length of all the 
y-axes. Each individual 
axis length will equal 
(YLEN/NZ) -.25 YLEN 
must be < 10. 
41-50] Constant CONVAR | Float~|Range value for all 
parameter ing curves. 
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NT I 


VARTABLE NAME 
OR CONTENTS 


| CARD | CARD 
TYPE | COL. 


51-80 | Not used 


5c =, 1-10 


| Number of 
different DELT's 


11-20 | Minimum 
| time value 
| 


| 21-80 | Not used 
5D 1-10 Number of times 
to use the DELT 
value found on 
this card 
11-20 | Increment 
| between time 
values. 
| 21-80 Not used 
SE 1-10 Number of 
curves 
' 11-80 | Not used 
1 SF 1-10 | Theta value for 
Ist curve 
11-20 | Theta value for 
2nd curve 


21-30 | Theta value for 
3rd curve 


PROGRAM 
NAME 


NDELTS 


TMIN 


NDELT 


DELT 


NZ 


TYPE 


Inte- 
ger 


Float- 
ing 


Inte~ 
ger 


Float- 
ing 


Inte- 
ger 


Float- 
ing 


Float- 
ing ¢ 


Float- 
ing 


DESCRIPTION 


Blank 


This controis the number 
of type 5D cards to be 
read 


First time value, basis 
for calculating all 
other time values. May 
be source or retarded 
time (see TCOL be'ow). 


Blank 

Each time value for the 
x-axis is equal to the 
previous value plus DELT. 
NDELT tells how many times 
to add DELT. 

See description of NDELT 
above. 

Blank 

Determines number of 

theta values to appear 

on card 5F below. Must 

be < 10 in present version. 


Blank 


Theta value in degrees 


Theta value in degrees 


Theta value in degrees 


CARD | VARIABLE NAME PROGRAM 
COL. | OR CONTENTS NAME 


31-40 | Theta value 


in degrees 
| for 4th curve 


| 
{ 
41-50 | Theta value | 2(5) value in degrees 
| for 5th curve | 
51-60! Theta value i 2(6) value in degrees 
for 6th curve | 
61-70 | Theta value | value in degrees 


for 7th curve 


71-80 | Theta value 


for 8th curve 


value in degrees 


EMP field to 
be plotted 


Must equal BPHI, ETHE, 
ERAD, or SIG7 


Not used Blank 


Time indicator TCOL=0, retarded time 
TCOL4¥0, source time 


x-axis spacing 
control and 
theta redefi- 
nition control 


INCX=1, spaced according 
to data; INCX=2, equally 
spaced, If INCX=3, a 
new set of 5E-5F cards 
are read to redefine the 
theta values. 


y-axis scaling ISCALE 


contro] 


ISCALE=0, scale all curves 
as one set; ISCALE4 0, scale 
each curve separately 


Blank | 


A standard deck for type 5 plots will consist of one type 5A, one type 5B, 


Not used 


one type 5C, as many type 5D as specified by the NDELTS parameter, one type 


5E, one type 5F (a second type SF is required if NZ is greater than eight), 


and as many type 5G cards as desired. If the last type 5G card sets 


INCX=3, then the series 5&-5G may be repeated. The last type 5G card must 
be followed by a blank card which signals the end of that set of type 5 
plots. 

Plots of types | through 5 may be requested in any order desired as 
long as all the cards needed for a particular type are grouped together. 
To stop execution of the plot cycle and close the plot file, a card with 


the word STOP in columns 1-4 must be placed after the last blank card. 
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I TEE a I 


VW. 


Sample Output 


Figure 
] 


2 


3 


Plot Type 
1 


1 


Descripcion of Sample Output 


Time Values 


Retarded 


Retarded 


Source 


Source 


Retarded 


Retarded 


Source 


Source 


Source 


Source 


Retarded 


Retarded 


Source 


Source 


Retarded 


Retarded 


Source 


Source 


Y-axis Sca'ing 


As one set 
As one set 
As one set 
Separately 
As one set 
As one set 
As one set 
Separately 
Separately 
Separate!v 
As one set 
Separately 
As one set 
Separately 
As one set 
Separately 
As one set 


Separately 


X-axis Spacing 


By value 


Equally spaced 


Equally spaced 


Equally spaced 


By value 


Equally spaced 


Equally spaced 


Equally spaced 


By 
By 
By 
By 
By 
By 
By 
By 
By 


By 


actual value* 
retarded value 
value 
value 
value 
value 
value 
value 
value 


value 


*Note: The actuai value data points should have the same relative spacing 
The following modification to program 


etna atte, etait netted ten aon tenet 
ae 


as the retarded value data points. 
COMPARE should eliminate the discrepancy: 
statement 1220 should read X(I) = I*RINC. 
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wee ene nee pee NEY RETNA MEE TON, BOER 


B-PHI 
(IN 1073 WEBER/Mam2) ALL curves are FoR R=7.SSe10! 
0.00 t=1.00010°° 
-0.76 . 
“1-52 
“2.27 
-3.03 
0.00 
-0.76 
-1.52 
-2.27 
-3.03 
0.00 
-0.76 
-1-S2 
-2.27 
-3.03 
0.00 
-0.76 
-1.52 
-2.27 
-3.03 
0.00 T=3.00"10-8 
-0.76 
-1.52 
-2.27 
~3.03 
0.00 
-0.76 
-1.52 
-2.27 
-3.03 
0.00 
-0.76 
-1.52 
-2.27 
-3.03 
0.00 eer caterer aren 
-0.76 
“1.52 
-2.27 
-3.03 
0.00 
-0.76 
-1.52 
-2.27 
-3.03 


0.00 Beh eG Top ooo 


.=).50010°8 


_., T22. 0081 0-8 


T=2.50"10°8 


T=3.50"10-8 


t=4.00010°8 


T24.50m10-8 


". 


t=5.000107% 


0.76 
71.52 
2.27 
-3.03 


U 


o 
Q 
o 
Q 


000° Zé 
oaga:29 
aaa‘ 2t 


aaa: 707 
goo° ag} 
o00°8?7 
Qoa° at 
gao-ae) 


Figure |. 
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B-PHI 


CIN 10° WEBER/MwsZ! ALL curves are FoR R=7.55810! 
0.00 T-1 0081078 
-0.76 “v2 
-1.52 


T=] 5010-8 


7=2.00"1075 


T=2.50"1078 


T=3.00010-8 


T=3.50"10°8 


T=4.00m1078 


tT=4.50m1078 


0.00 +... Saat ARE ERED nian ene Seemann Abaaes anne aamneen manne mente t=5.00"10°° 
-0.76 
-1.82 
-2 27 
-3.03 
OO a pene rane rf meen fon mare menpnreonrrefemer preempt T=S.50m10°° . 
-0.36 
~1.52 
2.27 & 
-3 03 
Ln JannEnnnT Sannin eet amas Sane Jeatiemets aelienent tiemariat nian Sennen toeeteei Renn Raenee eeeees Se ee a, 
o w n v ne] @ w aw ) a o ira) a oO ~ ~ ~ ~ am 
. NN nN Ld oO ~~ @ © wo Qo Q Q > un Q Q ce ~ wor 
<1 . . : : : = @ o @ om 
Q o oS i=] o ae) Q bad wo o Nn o ~~ ~ - . - . _™ 
Q Q io] i=] =] w un i=] i=] Q in] a un nH Qo o o Qa ov 
Q Q Q a a =] =) So o a Qo Q Q Q i=] Qa Qa i=] 
Q o a o o 
Figure 2. 
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7.5510! 
T=2 .5261077 
2.62"10-? 
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2.82010 
2.92810 
3.02810 
3.12"10 
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T=2. 
T=2. 
T= 
T=2. 
T=3. 
123. 
T=3. 
T=3. 
T= 


THETA 
(80.000 


£48 .000 


{18.000 


(08 .O0D0 


ALL CURVES ARE For R 


t0t .ODB 


95.750 


gt .750 


90 .600 
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Figure 3. 
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0.09 Ee es ee a ee eT , T=2.82"10-7 
-0.89 EXP=4 
-1.87 P 5 
at ee 
-3.83 
T=2.92m1077 
EXP=4 
T=3. a a 


ele eee ee eee , T=3.12"10°7 
: 7 EXP=4 


0.39 inst a 133281077 
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Figure 5. 
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Figure 9. 
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Figure 10. 
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T IN NSEC 


Figure il. 
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Figure 15, 
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Figure 16. 
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Figure 17. 
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B. BUDDHA 
1. Introduction 

Program BUDDHA was originally designed to process a maximum of two 
data tapes according to record number, to pilot right and left hand axes 
on a linear scale, and to overlay curves! The prevent version of BUDDHA 
has been expanded to process multiple tapes according to record name 
and/or record number and to plot right and left hand axes on a log scale. 

2. Input Revisions 

There are three additional data card types and card type Be has 
three new variables defined. 

Type C is the tape card, type D is the record card, and type € is 
the log card. A data set consists of a tape card followed by as many 
record cards as specified on the tape card. When all the record cards 
have been read, another tape card is read and the cycle continues until a 
blank tape card is read. An axis card is then read which is followed by 
a log card and as many curve cards as specified on the axis card. When al] 
the curve cards are read another axis card and a log card are read and the 
cycle continues until two blank cards are read which indicate all input 


data has been processed. 


1. Jones, D. L. and D. H. Stump, ELECTRA, ORESTES, AND SUPPORTING GRAPHIC 
DISPLAY CODES, U.S. Army Mobility Equipment Research and Development 
Center, Fort Belvoir, Virginia, February 1971, pp. 169-195. 


2. Ibid, p. 169. 
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21-80 


Not used 


Not used 


Not used 


VARIAB’.E NAME 
OR CONTENT 


PROGRAM 


NAME 


KARD (NS ,3) 


DATE 


| RUN 


1ALT 


1Z(NS) 


ITAPE 


NDATES 
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TYPE DESCRIPTION 


Same as report descrip- 
tion 


Data tape number (1, 
2,3, etc.) 


Same as report descrip- 
tion 


Blank 


Name of record to be 
plotted 


Run number 


Altitude Number 


Light intensity 


Blank 


Tape number 
(1, 2, 3, etc.) 


Number of record names 
to be processed for a 
particular tape. 


Blank 


D 
Not used 
i Not used 
2 E 1-4 
} 5-8 Not used 
aon 9-10 
d 11-16 Not used 
a 17-20 
4 : 21-25 


- CARD CARD VARIABLE NAME PROGRAM 
TYPE COL. OR CONTENT NAME 


TDATE(I, 


ITAPE) 


NRUNS 


MAXALT 


MAXYVAL 


[LOGO 
(INC) 


IPRT 


ILINE 


DECADE 
(INC, 1) 
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TYPE DESCRIPTION 


Alpha- | Record name on data 
numer~ | tape, where | runs 


ic from 1 to NDATES. 
Blank 

Inte- | Total number of runs 

ger for a particular 


record name 


Inte- | Total number of altitudes 


ger for a particular record 
name 
Inte- | Nuuber of maximum value 
ger records 
Blank 


Alpha~ | LOGO indicates a log 
numer- | axis 
ic 


Blank 


Alpha- | Print control for 
numer- | debugging 
ic 


Blank 
Alpha- |] Grid line indicator 
numer- 


ic 


Float- | Number of decades for 
ing left hand axis 


VARIABLE NAME 
OR CONTENT 


Not used 


PROGRAM TYPE DESCRIPTION 
NAME 


DECADE Number of decades for 


(INC, 2) right hand axis 


Blank 
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3. Data Revisions 

In the original program after all data cards were processed, the 
curve cards were ordered according to record number and the random access 
files were opened. The appropriate data tape was read and data was written 
on the corresponding random access file.? {in the present version, a tape 
card and as many record cards as specified on the tape card are read. For 
each record card that is read, the total number of records and the starting 
block number for that particular record name are calculated. The total 
number of records and the starting block number are then used to calculate 
the appropriate record number that corresponds to the record name on that 
curve card. The data tape number on that curve card is multiplied by ten 
thousand and is then added to the record number. The curve cards are then 
ordered according to record number and the random access files are opened. 

The reading of the datx tape or tapes has one modification. Because 
multiples of ten thousand based on the tape number are added to the record 
number, all the records from a particular tape are processed and then a new 
tape is processed. Data from odd number tapes are stored on one random 
access file and data from even number tapes are stored on Srotiet random 
access file. 

4, Plot Additions 
Program BUDDHA has been expanded to plot ordinate axis data on a 


log scale. The data points are plotted on a positive log axis, a negative 


3. Ibid., pp. 170-171 
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log axis, or on a split positive and negative log axis depending upon the 
behavior of the data. When the left hand axis array contains all positive 
data, the data points are scanned to determine the positive maximum and the 
positive minimum. If the number of cycles is defined by the user, the 
positive minimum is redefined as the difference between the positive max- 
imum and the number cycles. If the number of cycles is not defined by the 
user, the number of cycles is set equal to the difference between the pos~ 
itive maximum and the positive minimum. The negative maximum, the negative 
minimum, and the number of cycles are determined in a similar manner when 
the data array contains all negative data. Both sets of variables are 
determined when the array contains positive and negative data. The process 
is identical for the right hand axis data if appropriate. 

The length of the log axis is determined by the behavior of both 
the right and left hand axis data. The axes are 8.5 inches if the data 
arrays contain all positive or all negative data. If one array contains 
positive and negative data, the axes are split into a positive axis and a 
negative axis which are 4.18 inches long. The iog axes are drawn using 
subroutine LOGPLT (see 1.B.6). The time axis is located at the middle of 
the plot when the data is positive and negative and at the bottom when the 
data is either all positive or all negative. The time array is annotated 


with the minimum and maximum cutoff times us in the original program. * 


4, tbid., pp. 171-172 
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Ree 


& 


mers 


a ee 


in order to begin plotting, the data array is scanned to determine 


Oo ey eee 


the starting point and the end point of the positive and negative data for 


4 a i i 


F a particular curve. Any zerces that are encountered are omitted. If the 


data oscillates positive and negative more than ten timez, the variable 
is not plotted and message is printed. The logs of the data points are 
taken, converted into inches, and are plotted using subroutine LINES (see 
1,B.5) until all curves are drawn. Any points that are less than the 
defined minimum are set equal to the minimum. The ordinate axis is 


labeled and a title is written at the bottom of the plot. The right hand 


axis is then labeled if appropriate. An enclosing box is drawn around 


a . the plot and the paper is advanced until all plots are made, after which the 


| plot file is closed. 
aa : 5. Subroutine LINES 
. 4 Subroutine LINES is a plot routine designed for drawing data curves 


and annotating plots. Given two data arrays whict are to be plotted against 


Yai each other, subroutine LINES performs a curve fit to the data in a manner 
i . such that the curve can be represented by a series of equally spaced dots, 
4 line segments or both. Given the X-coordinate and the Y-coordinate, sub- 


routine LINES draws the appropriate symbol at the bottom of the plot. 


aii 
4 . 6. Subroutine LOGrLT 
§ 
a ; Subroutine LOGPLT is a special plot routine designed to draw and 
- i annotate positive and negative log axes. The log axis is annotated with 


tic marks at the power of ten and intermediate tic marks at 2, 4, 6, and 8 
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between the powers of ten. For negative data the power of ten is preceded 
by a negative sign. The tic marks are perpendicular to the axis line and 
are bisected by it. A positive angle, in degrees, produces a positive log 
axis and a negative angle produces a negative log axis. A switch, IFR, is 
used to determine a right or left hand axis. If IFR is not equal to one 
a left hand axis is drawn and if IFR is equal to one a right hand axis is 
drawn. 
7. CALCOMP Subroutines 

There are three subroutines from the CALCOMP plotting package that 
are called by subroutine LOGPLT. The subroutines are NUMBER, PLOT, and 
SYMBOL. 


NUMBER - writes floating point numbers. 


PLOT ~ converts pen movement from inches to actual plotter 
commands. 
SYMBOL - writes alphanumeric text and special symbols at any 


angle desired. 
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C. GREEK 
1. Introduction 
ce Subroutine GREEK is the controlling routine for a set of three 
subroutines which together draw the letters of the Greek alphabet using 
the CALCOMP plotting system. Subroutine GREEK can draw all 24 Greek 
letters in either upper or lower case. The user may specify any height 
and any orientation ang!e within the limitations of the plotting equip- 
ment. 
2. Calling Procedure 
Subroutine GREEK is called in Fortran by the following statement: 
CALL GREEK(XPOS, YPOS, SCALE, LETTER, THETA, ICAP) with the 
arguments defined as follows: 


XPOS - the X position (in inches) of the lower left-hand corner 
(before rotation) of the letter. 


YPOS ~ the y position (in inches) of the lower left-hand corner 
(before rotation) of the letter. 


SCALE - the scale height(in inches) desired for the letter. Some 
lower case letters will be drawn slightly smaller to main- 
tain proportionality between letters. A height which is a 
multiple of 31 times the plotter increment is recommended 
for best results. 


LETTER - an irteger number (1 to 24) specifying the position in the 
Greek alphabet of the letter to be drawn. 


THETA - orientation angle (in degrees) at which the letter is to be 
drawn. 0° - letter is vertical: positive angle ~- rotates 
letter counterclockwise. 


ICAP - selector for upper or lower case letters. ICAP =0 -- 
lower case, ICAP=] -- upper case. 
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Subroutine GREEK was written with the assumption that the user will 
open the plot file and establish an origin prior to calling the sub- 
routine. 
3. Data Control 
The data points for each letter were selected by first drawing 
the letter on a 31 x 31 grid and then approximating all curved lines by 
straight line segments. The coordinates of the starting and stopping 
points of each line segment were then recorded along with the pen 
up-down commands. Next, in order to minimize the amount of computer 
storage space needed for the subroutine, the data points were packed such 
that each 60 bit word contains the x and y coordinates and pen commands 
for five end points. When a letter is being drawn, the particular data 
points for that letter are then unpacked by subroutine SEPAR8 described 
below. 
4, Subroutine SEPAR8 
Subroutine SEPAR8 unpacks the 60 bit computer words and separates 
the data into x and y coordinates and pen commands. [t performs this 
operation by using a series of left and right shifts and masking functions. 
5. Function Subroutine |SHIFT 
Function subroutine ISHIFT performs left or right shifts of | to 
60 bits. It is written in COMPASS (Control Data Corporation assembly 


language) for use on a CDC 6000 series computer. 
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6. Special Considerations 
. Since the data points have been packed intc 60 bit words and 
function subroutine ISHIFT is written in COMPASS, the subroutine GREEK 
package presented here can only be utilized on a CDC 6000 series computer. 
The package requires 633 (octal) words of core and utilizes the CALCOMP 


routine PLOT. 
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7. Sample Output 


(height of letters = 0.31) 


118 


Nea 


rAd 


Ce hoe Ge GR ee 


Vive OMG eS Be 26 VO 


AB iE. Onl Ra 


Ne Ose 2. Oe a 
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It. MENCR MODIFICATIONS 
A. ELECTRA 
1. Description of Modifications 

Program ELECTRA H52 the ground burst EMP code, received two 
minor modifications since it was last documented. All input/output (1/0) 
operations using tape or drum were revised in order to cut running time 
and to make more efficient use of data tapes. 

The running time was cut by replacing all lengthy 1/0 lists, which 
used implied DO loops, with a variably dimensionea subroutine (IOSAYE) 
which can write or read any length array as a single operation. This 
modification should produce up to 90% reduction in running time for the 
operations which it affects. 

More efficient use of data tapes was accomplished by revising the 
contents of the data output array (OUTARR). By removing some unnecessary 
data items from the array and repositioning other items, the dimensions 
of the array were reduced from 2ix7 to 21x4 resulting in a decrease of 63 
words. The words in the reduced array were then packed two into one by an 


assembly language packing routine (PACK). The net result was to rediice 


1. Jones, D. L. and D. H. Stump, ELECTRA, ORESTES, and Supporting Graphic 
Display Codes, U.S. Army Mobility Equipment Research and Development 
Center, Fort Belvoir, Virginia, February 1971, pp. 1-36. 


2. Borbely, J. A. and D. L. Jones, ELECTRA, An Electromagnetic Pulse 


Fortran Program (User's Guide), U.S. Army Mobility Equipment Research 
and Development Center, Fort Belvoir, Virginia, October 1969. 
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the numoer of words written to tape, at a given time step and range step, 
from 147 to 43. Uf the number of angles used for a given run is less 


then 18, the number of words written to tape will be reduced even more. 
A 
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B. ORESTES 
1. Description of Modifications 

Program orEsTES! received two minor modifications since it was 
last documented. All input/output (1/0) operations using tape or drum 
were revised in order to cut running time and to make more efficient use 
of data tapes. 

The running time was cut by replacing all lengthy I/0 lists, 
which used implied DO loops, with a variably dimensioned subroutine 
(1OSAVE) which can write or read any length array as a single operation. 
This modification should produce up to 90% reduction in running time for 
the operations which it affects. 

More efficient use of data tapes was accomplished by revising the 
contents of the data output array (OQUTARR). By removing some unnecessary 
data items from the array and repositioning other items, the dimensions of 
the array were reduced from 27x11 to 25x9 resulting in a decrease of 72 
words. The words in the reduced array were then packed two into one by an 
assembly language packing routine (PACK). The net result was to reduce 
the number of words written to tape, at a given time step and range step, 
from 297 to 114. If the number of angles used for a given run is less than 


24, the number of words written to tape will be reduced even more. 


]. Jones, D. L. and D. H. Stump, ELECTRA, ORESTES, and Supporting Graphic 
Display Codes, U.S. Army Mobility Equipment Research and Development 
Center, Fort Belvoir, Virginia, February 1971, pp. 37-69. 
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C. REDACT 
1. Description of Modifications 

Program REDACT has been modified and is now a subroutine with a 
driver program DATAPK. This modification was made in order to process 
data tapes from the modified ELECTRA which have a packed OUTARR (see II. 
A.1.). 

DATAPK opens and closes the plot files and calculates the number 
of words that are packed in CUTARR. In subroutine REDACT, the column and 
row indicatois for data editing have been modified to correspond to the 
modifications made in ELECTRA. In order to process the data tape, the 
column and row indicators for the particular variable to be plotted are 
used to determine the word to be unpacked. The data array is unpacked by 


an assembly language routine (UNPACK) . 
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D. REDACTO 
t. Description of Modifications 

Program REDACTO has been modified and is now a subroutine with a 
driver program DATAPKO. This modification was made in order to process 
data tapes from the modified ORESTES which have a packed OUTARR {see 11.B.1). 

DATAPKO opens and closes the plot files and calculates the number 
of words thac are packed in OUTARR. In subroutine REDACTO, the cotumn and 
row indicators for data editing have been modified to correspond to the 
modifications made in ORESTES. In order to process the data tare, the 
column and row indicators for the particular variable to be plotted are used 
to determine the word to be unpacked. The data array is unpacked by an 
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assembiy language routine (UNPACK). 
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